Coupling of the lattice and superlattice deformations and hysteresis in thermal 
expansion for the quasi one-dimensional conductor TaS3 
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An original interferometer-based setup for measurements of length of needle-like samples is de- 
veloped, and thermal expansion of o-TaS3 crystals is studied. Below the Peierls transition the 
temperature hysteresis of length L is observed, the width of the hysteresis loop SL/L being up to 
5 ■ 10~ 5 . The behavior of the loop is anomalous: the length changes so that it is in front of its 
equilibrium value. The hysteresis loop couples with that of conductivity. The sign and the value of 
the length hysteresis are consistent with the strain dependence of the charge-density waves (CDW) 
wave vector. With lowering temperature down to 100 K the CDW elastic modulus grows achieving 
a value comparable with the lattice Young modulus. Our results could be helpful in consideration 
of different systems with intrinsic superstructures. 

PACS Numbers: 71.45. Lr, 65.40. De 



Internal degrees of freedom is a feature of a random 
system; in principle, they can give rise to metastable size. 
states resulting, say, to hysteresis in thermal expansion.El 
A special class form the compounds with intrinsic su- 
perstructures. Comprising two periodicities, gener- 
ally incommensurate, the compounds occupy intermedi- 
ate place, between genuine aperiodic and truly periodic 
systemsE In these systems, such as charge- and spin- 
density waves (CDW and SDWi,H Wigner cristals, super- 
conductors in magnetic fieldsjj structurally incommen- 
surate crystal phases,u the superstructure periodicities 
could be varied by external fields or temperature changes. 
The resulting metastable configurations can be reflected 
back oato-the elastic properties and size of the underlying 
lattice,ETEl though this question is still poorly understood. 

Quasi 1-dimensional conductors with CDW belong to a 
widely studied class of materials, in which intrinsic super- 
structure develops through the Peierls transition^. When 
electrons condense into CDW they form a deformable 
medium, - an electronic crystal. Deformation of the 
CDW affects their main static and dynamic properties 
and gives rise to metastability and hysteresis. 

The straightforward treatment of the CDW as a spring, 
whose strain is just applied to the crystal at the ends or 
via the impurities is not valid. Moreover, in the simple 
one-dimensional model the strains of the CDW and the 
crystal do not couple at all: if initially the CDW are re- 
laxed, any change of the crystal length would not draw 
the CDW away frorruthe equilibrium, i.e. give rise to 
a CDW deformation,0 as it was noticed in Refs. |^-|TT|. 
Similarly, once the CDW is deformed, any change of the 
lattice constant, c, would neither decrease nor increase 
the deviation of the CDW wavelength A from the equi- 
librium value, X eq . So, within this model a CDW defor- 
mation would not give rise to a length change. 

At the same time, the interaction of the CDW and 
the lattice is clearly seen from the elastic anomalies_ipi 
eluding a drop p£ the Young modulus of the latticeJjtiJ 
Yi, up to 4%Ji3 when the CDW become depinned. 



MozurkewichB has concluded that the lattice deforma- 
tion does give rise to a deformation of the CDW: when 
the CDW are at rest, they cannot relax, and so they 
contribute to the total elastic energy; the sliding CDW 
relax rapidly and their contribution drops outJ3 For his 
model Mozurkewichu introduced an empirical parameter 
g reflecting the deviation of \ eq from the simple one- 
dimensional model: 



SX eq (c)/X eq (c eq ) = (g+ l)dc/c e 



(1) 



g 0. However, no way was proposed to estimate g. 

Another remarkable observation was reported in 
Rcf. p2t the length L of TaS3 samples as a function of 
electric field demonstrated hysteresis partly scaling with 
that of resistance. The values of L obtained at differ- 
ent directions of the voltage sweep differed by ~ 10~ 6 . 
This result was also treated in terms of coupling of the 
CDW strains with the deformation of the pristine lattice. 
Note that the field-induced length change (as well as the 
change of resistance) results only from an inhomogene- 
ity of the sample properties: obviously, the length could 
depend on the voltage polarity osly if the inversion sym- 
metry of the sample is brokentd Thus, study of the 
electric-field induced hysteresis of length can not provide 
complete understanding of the CDW-lattice coupling. 

To study the effect of the CDW deformation on the 
crystal length it could be more fruitful to observe thermal 
hysteresis of L. With changing temperature the g-vector 
(q = 2tt/X) falls behincLits equilibrium value, which is 
temperature dependent ;li3 the deviation of q from the 
equilibrium is limited by the critical deformation at which 
the phase slippage (PS) begins.E-3 So, thermal cycling 
creates a CDW deformation, which is relatively uniform 
along the sample and is of the maximum possible value. 

For our study we chose o-TaS3 as a representative 
CDW conductor. Among the quasi one-dimensional con- 
ductors TaS3 is one of the most widely studied, includ- 
ing the elastic properties. TaS3 demonstrates the Peierls 
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transition at Tp — 220 K, below which the resistance fol- 
lows an activation law revealing the half-gap A = 700 K. 
The dependence R(T) demonstrates a pronounced hys- 
teresis, the width of the loop 5R/R being up to 50 % at 
T around 100 K. 

In the present Letter we report temperature hystere- 
sis of length for o-TaS3; the hysteresis couples with that 
of resistance having an anomalous sign. A quantitative 
treatment of the effect is proposed. 
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FIG. 1. a) Arrangement of a TaS3 sample on the substrate 
and the scheme of the reflections and interference of the laser 
beams (the incidence is close to normal), b) An example of 
an image obtained, c) - the section of the image, d) - the 
resulting profile of the sample. The open circles correspond 
to the minima in the section (Fig.lc), and the closed circles 
— to the maxima. 

We studied samples of TaS3 with typical length 1 mm, 
width 15-35 /im, and thickness 5-10 times less. The 
samples were arranged on a transparent glass substrate 
(Fig. la). The contacts were fixed on the substrate with 
indium, while the central part of the sample formed an 
arc. The inner surface of the substrate played the role of 
a semi-transparent mirror. The laser beam with wave- 
length A; = 635 nm fell down through the substrate, and 
partially reflected from the inner surface of the substrate 
and from the sample surface forming an interference pat- 
tern (similar with the wedge interference). The pattern 
was fixed with a video camera combined with a micro- 
scope. An example of such a pattern is shown in Fig. lb. 
The neighboring dark or bright fringes correspond with 
change of the sample distance y from the substrate sur- 
face by ±A;/2. Making sections of such images (an exam- 
ple is shown in Fig. lc) we obtained the profiles y(x) of 
the samples, like the curve presented in Fig. Id. Finally, 
the length of such curves was calculated. We detect a 
length change about 5 A which corresponds to the rela- 
tive length change 5L/L — 5 • 10~ 7 , though in principle 



the sensitivity could be improved.0 Our estimates have 
shown that the contributions to 8L due to the changes of 
stress and Yj are negligible Jl3 Thus, the observed changes 
of L with temperature were associated with the change 
of the equilibrium sample length with respect to the sub- 
strate. To have absolute results we have measured ther- 
mal expansion of the glass substrate in addition. 
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FIG. 2. a) Temperature dependence of SL/L with respect 
to glass. The absolute thermal expansion is shown in the 
inset, b) R(T) measured simultaneously, c) The width of the 
length (the circles) and conductivity (the solid line) hysteresis 
loops as a function of T. d) The ratio (5L/ L)/(8X/X) (open 
circles) and gY c /Yi (black circles) resulting from Fig. 2c. We 
use Eq. ^ and the fi(T) dependence.!^ The solid line indicates 
Y c oc exp(470/T). 

Figs. 2a) and b) show the temperature dependences 
of length (with respect to glass) and of the resistance 
for the representative sample of TaS3. In the inset to. 
Fig. 2a) the absolute dependence of length is shown.EZl 
The dependences of length and resistance both clearly 
show hysteresis below the Peierls transition temperature 
Tp = 214 K. The length hysteresis loop opens imme- 
diately below Tp. The development of hysteresis is the 
dominating effect near Tp , on whose background it is dif- 
ficult ±o distinguish a feature coupled to the transition in 
itself.ta The length hysteresis grows below Tp achieving 



2 



5L/L = 5 ■ 1CP 5 (Fig. 2a), which by 1.5 orders of magni- 
tude exceeds the maximum value reported in Rcf. [l2|. At 
the same temperature the states approached from lower 
temperature have higher length than those approached 
from above. Correspondingly, application of electric field 
exceeding the threshold reduces the length after heating 
the sample, and increases after cooling. Thus, the main 
effect of electric field is the relaxation of the thermally 
induced states. The length change induced by the field 
in itself was ~ 10 -6 , like in Ref. [l^. Note, that the sign 
of the hysteresis is non typical: the length is in front of 
its equilibrium value (inset to Fig. 2a), but not behind. 

The length hysteresis apparently couples with the 
resistance hysteresis, higher length corresponding to 
higher resistance at a given temperature.Ej The cou- 
pling between R and L is well seen from Fig. 2c, where 
the temperarure dependence of SL/L is presented to- 
gether with the loop of conductivity, 5a = l/R CO oi — 
1/Rheat- Note that iu the unipolar approximation 5a ks 
03Oo((5A/A)(/z/^3oo),E3 where the index "300" marks the 
room-temperature values. So, the a hysteresis reflects 
that of the CDW wavelength to the accuracy of the tem- 
perature dependence of the mobility fi, which in turn 
can be found from Ref. Fig. 2d shows the ratio 

(5L/L)/{5\/\) vs. T. At the low temperatures 5L/L 
is ~ 6 % of the CDW deformation. 

As we have noticed in the beginning, the interactions 
of the CDW and the sample cannot be presented as that 
of two springs connected in parallel (the lower inset to 
Fig. 3) . Though we do not see a way to obtain quantita- 
tively the strain dependence of the q- vector for TaS3, one 
can recall experimental results, from which it is possible 
to estimate the value of g (Eq. |l|). It has been foundE3 
that at a certain value of uniaxial strain S* the proper- 
ties of the CDW change abruptly: e.g., the non-linear 
conduction nearly disappears. This has been attributed 
to the transition of the CDW to the 4-fold commensu- 
rability. —Comparing the normalized dependences of the 
g-vectoiO and of the strain S*t3 on tecnperature (Fig. 3) 
one can see that they become similarE3 if one multiplies 
the value of S* by 6. Thus, one can assume that a defor- 
mation of the lattice induces the change of the equilib- 
rium CDW wavelength X eq (c) in accordance with Eq. (|l|) 
with j«6. Such a large jaalue of g indicates that the sim- 
ple 1-dimensional modelEI cannot even roughly describe 
the strain- induced change of X eq . Evidently, the change 
of X eq is dominated by transverse effects: the longitudi- 
nal strain decreases the thickness of the sample resulting 
in an increase of the intrachain coupling. This modifies 
the form-,of the Fermi surfaces, and consequently, the q 
-vector.E3 To comprehend the sign and the value of the 
q change one should study the evolution of the Fermi 
surface in detail. In principle, we cannot exclude that 
g is temperature dependent. However, as we shall see 
below, the assumption that g = const, namely g w 6, is 
consistent with our experiment. 

For the next step, we present the elastic energy density 
W as a sum of the lattice and the CDW energies: 
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FIGj-3. The temperature dependences of S* — S(0) (open 
circles)Ei, and of the wave vectoitj (q-q(0))/q(0)) (closed 
circles) .Ej S* is multiplied by g = 6. Insets: the equivalent 
schemes illustrating the interaction of the CDW (the thin-line 
spring) and the pristine lattice (the thick-line spring). The 
lower-left sketch is the naive scheme, which is invalid. The 
upper-right sketch illustrates Eq. [| (g > 1). 

w = krtC—^) 2 + r c ( A ~ Aeg(c) ) 2 ), (2) 

where Y) and Y c are the elastic moduli of the lattice and 
the-CDW respectively. Taking into account the condition 
(|l])c-3 and assuming that no PS occurs (A/c = const), we 
can minimize W and obtain the resulting length change: 

5L_5c__ Y c 5X Yc5X 
L ~ c eq ~ 9 Y l +g 2 Y c A ~ 9 Y x A' (3) 

where 5\ is the initial CDW deformation (at fixed c). 
The approximation implies that g 2 Y c <C Y\. Eq. ^|is quite 
transparent: the crystal deforms as a spring connected in 
parallel to the CDW, but with a factor —g. g > means 
that, say, compressed CDW (A < \ eq ) would result in 
a decrease of the sample length, which agrees with our 
observation. E.g., cooling corresponds to growth of A e<3 lii2, 
so the CDW is in compressed state (A — \ eq < 0), which 
results in decreased L (Fig 2a). The CDW - crystal 
interaction could be illustrated with a scheme consisting 
of two springs connected via a lever (Fig. 3, upper inset). 

From Eq. | it follows that the ratio (5L/L)/(5X/X) 
(Fig. 2d) equals gY c jY\. Note, however, that with g = 6, 
g 2 Y c at low temperatures is comparable with Y;, and one 
should use the exact version of Eq. |[ The resulting value 
of gY c /Yi is also shown in Fig. 2d. One can see the growth 
of the CDW elastic modulus with loweidag temperature, 
though slower than Y c oc exp(A/T)EJE3 (see the solid 
line in Fig. 2d).C3 Note that the curve (5L/L)/(5X/\) vs. 
T resembles the temperjaiure dependences of the lattice 
softening in electric ficldlij. This could be expected, since 
from Eqs. |l| and ||it follows that the value g 2 Y c is added 
to Yi when the CDW are at rest, while in the sliding state 
the CDW contribution drops outcl. 
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From the condition of neutrality!^! accompanying the 
CDW deformation we can obtain for the CDW Young 
modulus: Y c = ^^0, where ^ is the derivative of the 
chemical potential level by the CDW g-vector, and s is 
the area per chain. Substituting this in Eq. pi we obtain: 



5L 



25( 



L ~ g ~\sYi 



(4) 



For the unipolar (p-type for TaSs) conductiorJU'LH S( 
Tln(R heat /R coo i). Substituting into Eq. g — 6, 5C, = 
50 K (for T = 120 K, Fig. 2b), A = 10 X s = 80 A 2 
(Ref. |7|), Yt — 380 GPaEl, we obtain SL/L = 3- 10~ 5 , 
in agreement with the experiment (Fig. 2b). This agree- 
ment supports the approach we used to estimate g, and 
proves that under a certain uniaxial strain (S^Lthe CDW 
in fact transit to the 4-fold commensurability.E 2 ] 

In conclusion, we have observed thermal hysteresis of 
length for the quasi one-dimensional conductor TaS3. 
In the framework of the model proposed the observed 
anomalous sign and value of the hysteresis loop are con- 
sistent with the dependence of the q- vector on the lattice 
strain. The results could be helpful in understanding the 
behavior of other systems with intrinsic superstructures, 
whose parameters could be rearranged due to the tem- 
perature effects or external fields. 
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17301, 00-02-22000 CNRS), the State programs "Physics 
of Solid-State nanostructures" and "Low-Dimensional 
quantum structures" , NWO and INTAS-0 1-0474. 



1 For a recent example see: P. Nagel, V. Plaster, C. Meingast 
A. Rykov, and S. Tajima, Phys. Rev. Lett 85, 2376 (2000) 

2 J. Voit, L. Perfetti, F. Zwick, H. Berger, G. Margaritondo 
G. Griiner, H. Hochst, M. Grioni, Science 290, 501 (2000) 

3 J.W. Brill in: Handbook of Elastic Properties of Solids 
Liquids, and Gases, Volume II, Ed. by M. Levy, H.E. Bass 
and R.R. Stern, (Academic Press 2001), pp. 143-162. 

4 V. V. Eremenko, V. A. Sirenko, H. Szymczak, and A. Nabi- 
alek, Low Temp. Phys. 25, 225 (1999). 

5 H. Z. Cummins, Physics Reports 185, 211 (1990). 

6 Charge Density Waves in Solids, edited by L. Gor'kov, 
G. Griiner (Elsevier Science, Amsterdam, 1989); G. Griiner 
in: Density Waves in Solids (Addison- Wesley Reading, 
Massachusetts, 1994); 

7 This follows from the fact that on longitudinal deformation 
of a sample the number of electrons per the lattice period 
does not change, and A eq (c) changes proportionally to c, 



SX eq (c)/X eq (c 



5c/ c 



8 G. Mozurkewich, Phys. Rev. B 42, 11183 (1990). 

9 J.W. Brill and W. Roark, Phys. Rev. Lett 53, 846 (1988); 
Phys. Rev. B 36, 2969 (1987). 



10 R.L. Jacobsen and G. Mozurkewich, Phys. Rev. B 42, 2778 
(1990). 

11 Z.G Xu and J.W. Brill, Phys. Rev. B 45, 3953 (1992). 

12 S. Hoen, B. Burk, A. Zettl, and M. Inui, Phys. Rev. B 46, 
1874 (1991). 

13 Z.Z. Wang, H. Salva, P. Monceau, M. Renard, C. Roucau, 
R. Ayroles, F. Levy, L. Guemas, and A. Meerschaut, J. 
Phys. Lett. France 44, L315 (1983). 

14 V.Ya. Pokrovskii and S.V. Zaitsev-Zotov, Synthetic Metals 
32 321 (1989). 

The interferometric methods for thermal expansion resolve 
10" 2 - 5 • 10" 4 of a fringe [T.H.K. Barron, J.G Collins and 
GK. White, Adv.Phys. 29, 609 (1980)]; in our configura- 



tion this corresponds to SL/L ~ 10 



10" 



A.V. Golovnya, V.Ya. Pokrovskii, to be published. 

17 D. Maclean and M.H. Jericho, Phys. Rev. B 47, 16169 
(1993): the presented thermal expansion (at considerable 
sample strain and at high electric fields) is in rough agree- 
ment with our data. 

18 If one takes L(T) = (L cool (T) + L heat (T))/2, he will obtain 
a decrease of length ~ 10~ associated with the Peierls 
transition, in qualitative agreement with Ref. [d]. 

19 In Ref. Il3 the state with smaller L corresponded to Ipgher 
R. Evidently, both the R and L changes reportedE3 re- 
sult from non-complete compensation of the contact CDW 
deformations of opposite signs. In this case one can obtain 
an arbitrary relation between the R and L hysteresis domi- 
nated by the non- linear effects. Concerning our experiment, 
the small difference in the forms of the loops of L and R 
(Fig. 2c), as well as the small loop of L above Tp (Fig. 2a), 
is evidently within our experimental error. 

20 D.V. Borodin, S.V. Zaitsev-Zotov and F.Ya. Nad', Zh. 
Eksp. Teor. Fiz. 93, 1394 (1987); [Sov Phys. JETP 66, 
793 (1987)]; this relation could be easily obtained from the 
neutrality condition if one notes that a variation of the q- 
vector changes the number of the states under the Peierls 
gap, and, thus, the number of quasiparticles. 

21 Yu.I. Latyshev, Ya.S. Savitskaya, and V.V. Frolov, JETP 
Lett. 38, 541 (1983). 

22 V.B. Preobrazhensky, A.N. Taldenkov, I.Yu. Kalnova, 
Pisma v ZETF, 40, 182 (1984) [JETP. Lett 40, 944 (1984)]; 
V.B. Preobrazhensky, A.N. Taldenkov, S.Yu. Shabanov, 
Solid State Commun. 54, 1399 (1985). 

23 We took S(0) = 3.9 • 10~ 3 , and q(0) = 0.2503c* suggesting 
that both dependences follow the Arrhenius law (the slope 
of the line corresponds to 700 K, which is close to A). 

24 Eq. H implies that the principal part of the CDW elastic 
energy is associated with their longitudinal deformations, 
while Eq. (Q) takes into account the dependence of X eq on 
all components of the crystal deformation. 

25 S.N. Artemenko and A.F. Volkov, Zh.Eksp.Teor.Fiz. 81, 
1872 (1981);[Sov. Phys.JETP, 54, 992 (1981)] 

26 Actually the CDW deformation is spatially non-uniform, 
so one can expect non-linear effects, e.g. the reduction of 
Y c at low T. 

27 P.Monceau in: Electronic Properties of Inorganic Quasi- 
one-dimensional Conductors, Part 2. Ed. by P. Monceau. 
Dortrecht: D.Reidel Publ. Comp., 1985. 

28 J.W. Brill, Solid State Commun 41, 925 (1982). 



4 



